On generation of Crab giant pulses 
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Abstract 

We propose that Crab giant pulses are generated on closed magnetic 
field lines near the light cylinder via anomalous cyclotron resonance on 
the ordinary mode. Waves are generated in a set of fine, unequally spaced, 
narrow emission bands at frequencies much lower than a local cyclotron 
frequency. Location of emission bands is fitted to spectral structures seen 
by Eilek et al. (2006). 

To reproduce the data, the required density of plasma in the giant 
pulses emission region is much higher, by a factor ~ 3 x 10 , than the 
minimal Goldreich- Julian density. Emission is generated by a population 
of highly energetic particles with radiation-limited Lorentz factors 7 ~ 
7 x 10 7 , produced during occasional reconnection close to the Y point, 
where the last closed field lines approach the light cylinder. 
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1 Introduction 

Giant pulses (GPs) are relatively rare durations of intense radio outbursts that 
are clearly a special form of pulsar radio emission (e.g. Kuzmin 2007). Overall, 
the properties of GPs are (i) peak fluxes densities can exceed thousands (and, in 
case of Crab, nearly a million) of times the peak flux density of regular pulses 
(the Crab pulsar was discovered by observations of its giant pulses ( Staelin & 



Rcifcnstein 



1968 1); on average, energy coming in GPs can exceed the energy 
in an average pulse. Distribution of peak fluxes is power-law (unlike regular 
radio pulses which have Gaussian distribution); (ii) GPs typically have short 
duration from several microseconds down to few nanoseconds ( |Hankins et al 



2003[ ) (much shorter that sub-pulses in the average pulse) and emit narrowband 



radiation (Popov & Stappers 20031; at lower time resolution, overlap of many 
narrowband emission spectra may mimic a broadband spectrum ; (iii) GPs seem 
to be limited to the edges of main profiles (e.g. Popov et al. 2006[ ) (see, though, 
unconfirmed observations of GP in HFCs by Jessner et al. (2005)); (iv) GPs are 
seen in only 11 pulsars; (v) there is no correlation between GPs and emission 
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at higher frequencies (nothing in X-rays ( Lundgren et al. 1995 ) and a marginal 
3% correlation in the optical ( Shearer et al.[ 2003)); (vi) there seems to be no 
strong correlation between pulsar properties (period, magnetic field , magnetic 
field at the light cylinder, luminosity) and GP phenomenon. 

In addition, recent observation of Crab pulsar ( |Eilek fc Hankins |2007 1 iden- 
tified unique features of GPs associated with the interpulsc (IGPs): (i) IGPs 
spectra consist of a number of relatively narrow frequency bands; (ii) spacing 
between the bands is proportional, Av/v ~ 0.06; (iii) emission at different bands 
start nearly simultaneous, perhaps with a small delay at lower frequencies; (iv) 
sometimes there is a slight drift up in frequencies; all bands drift together, keep- 
ing the separation nearly constant; (v) emission bands are located at 4—10 GHz, 
continuing, perhaps, to higher frequencies, but not to lower frequencies; (vi) all 
IGPs show band structure. These are very specific properties that allowed us 
to build a quantitative model of pulsar radio emission (applicable, at least, for 
Crab GPs associated with interpulse). 



et al. 19851 



Theoretically, GPs were proposed to be nonlinear solitons ( Mikhailovskii 
Petrova ( |2004 1 proposed generation by induced scattering of 



waves leading to a redistribution of the radio emission in frequency. Weather- 



all (20011; Hankins et al. (20031 proposed strong plasma turbulence; Eilek & 



Hankins (20071 discuss a possibility of interference fringes. In spirit, our work 
is closest to Istomin (2004), who proposed that GP are generated on Alfven 
waves during a reconnection event close to the light cylinder. On the other 
hand, discovery of the banded structure of IGPs is inconsistent with emission 
on Alfven waves (see Appendix [AJ, and, in addition, Alfven cannot leave mag- 
netosphere and need to be converted into escaping modes. 

Before laying out the details we give here a brief overview of the model, 
Fig. [l] The model requires that a high density plasma, with density ~ 10 
times Goldrech- Julian density, is present on closed field lines. In addition, an 
occational high energy beam, with radiation-limited Lorentz factor ~ 10 , is 
propagating along the field. These requirement come from actual fitting of ob- 
served emission bands to a particular cyclotron resonance conditions. The fit 
requires that the bulk plasma be not moving relativistically (hence association 
with closed field lines), the value of magnetic field be the minimum possible in 
the magnetosphere (even the factor of two difference of magnetic field strength 
between the magnetic pole and equator are important), emission region be lim- 
ited to small volume (hence association with the Y-point). 



2 Wave dispersion in pulsar magnetosphere 



To guide the reader through the following derivations, we first re-derive the basic 
properties of wave dispersion in pair plasma in a strong magnetic field of pulsar 
magnetosphere. In the standard model of pulsar magnetospheres (Goldreich & 



Julian 1969 1 rotating, strongly magnetized neutron stars induce strong electric 
fields that pull the charges from their surfaces. Inside the closed field lines 
of the neutron star magnetosphere, a steady charge distribution is established, 
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Figure 1: Generation region of giant pulses in Crab. High density plasma is 
trapped on closed field lines near the light cylinder. Occasional reconnection jets 
produce high Lorentz factor beams that propagate along magnetic field lines and 
emit coherent cyclotron-Cherenkov radiation at anomalous Dopplcr resonance. 



compensating the induced electric field. It is typically assumed that on closed 
field lines, particle density is at a minimum and equals charge density, the 
Goldrcich- Julian density uqj = £IB /(2irce). In fact, the real particle density 
may and does exceed this minimum value, as we argue in this paper. There 
is already a clear evidence that in case of the double pulsar PSR J0737— 3039 
plasma density on closed field line exceeds the Goldreich-Julian density by a 
factor - 10 4 - 10 5 ( Lyutikov fc Thompson 2005 1, though it was not clear if this 
is specific to the double pulsar. One of the implications of the present model 
is that high density plasma is present on closed field lines of isolated pulsars as 
well. 

As a simple parametrization, we normalize the plasma density to the Goldreich- 
Julian density, n p = A uqj. (We neglect the thermal motion of plasma particles 
and drifts resulting form curvature of field lines and inhomogenuity of magnetic 
field). The surface polar magnetic field of Crab pulsar is 3.7 x 10 12 G. Below we 
will argue that Crab CPs are generated near a magnetic equator, where mag- 
netic field at a given radius is two times smaller. Since the resonant frequency is 
proportional nearly to the third power of the local magnetic field ([8| , this factor 
of 2 is important. Thus, near the light cylinder the cyclotron frequency at the 
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equator is 

9-5 x 10 12 _! 

For the density parameter A = 3.3 x 10 5 (see below), the corresponding plasma 
frequency is 

3 4 x 10 10 

uj p = \/2\£Ilob = — r^rads -1 ( 2 ) 
(r/R LC y /2 

{y p = ujp/{2ir) = 5.5 x 10 9 Hz). Thus, near the light cylinder uj p /ujb ~2x 10~ 3 
(this is the maximum value of this parameter as a function of r. 

For e ± plasma in magnetic field the dispersion relation factorizes giving two 
modes: the X mode with the electric vector perpendicular to the k-B plane 
and two branches of the longitudinal-transverse mode, which we will call L-0 



and Alfven waves, with the electric vector in the k-B plane (Arons & Barnard 
1986 see Fig. [2j. X waves is a subluminal transverse electromagnetic wave with 



a dispersion relation 



2 



n 



2 = i-^K (3) 

here n = kc/u is refractive index, u>b = eB/mc is cyclotron frequency, lo v = 
\J 'Ami±e 2 1 'm is a plasma frequency of each species (so that for pair plasma the 
total plasma frequency is V2u> p ) . The Alfven -L-0 mode satisfies the dispersion 
relation 



,2 



(uj 2 - 2ujT){^ - ui%) - 2co 2 -- 2 



tot sin 



(4) 



Alfven branch is always subluminal while L-0 mode is superluminal at small 
wave vectors and subluminal at large wave vectors. 

2.1 L-O mode 

As we argue that the GPs are generated on the L-0 mode, let us discuss its 
properties in more detail. The L-0 mode exists in a frequency range 



V2uj p < uj < (uj% + 2uj p + (uj 2 b + 2uj p ) 2 — 8lu b oj 2 cos 2 ~\ 



^ 2 ^ 2 sin 2 i 

KUOB+— 

It has a reflection point at y2u> p and resonance at the upper bound. L-0 mode 
becomes luminal at frequency ujq = \j2ui 2 + uj 2 b sin 2 9. Its polarization is 

e x _ {uj 2 -uj 2 b ){uj 2 -2uj 2 p )cote 



e z uj 2 (uj 2 -u% -2lo 2 ) 



(6) 



(for wave vector in x — z plane), so that the L-0 mode is nearly longitudinal 
for fee <C too (this becomes a Langmuir wave for parallel propagation) and is 
nearly transverse for kc 3> luq. At the luminal point ujq polarization is quasi- 
longitudinal for 8 -c (ujpjujB) 2 and is quasi-transverse for larger angles. 
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3 Generation of giant pulses with banded fre- 
quency structure 

3.1 Anomalous cyclotron resonance on L-O mode 

The anomalous cyclotron resonance condition is 

u) — knvti = su>b/j, for s < (7) 

where k\\ and v\\ are components of wave vector and particle's velocity along 
magnetic field , 7 is Lorentz factor of fast particles and s is an integer. It is clear 
that the necessary requirement for the anomalous cyclotron resonance is that 
the refractive index of the mode be larger than unity and that the parallel speed 
of the particle be larger than the wave's phase speed. The physics of emission is 
similar to the Cherenkov process, except that during photon emission a particle 



increases its gyrational motion and goes up in Landau levels (Ginzburg 19851. 
The energy is supplied by parallel motion. Importance of anomalous cyclotron 
resonance for pulsar radio emission has been discussed by |Machabeli !k Usov| 



(19791; Kazbegi et al. (19911; Lyutikov et al. (19991. In these papers, emission 



was argued to be generated on open field lines at large distances from the start. 

The resonant frequencies corresponding to anomalous cyclotron resonance 
for the X and L-0 modes for our fiducial parameters (sec below) arc shown in 
Fig. [3j Overall, L-0 mode dispersion is quite complicated, depending sensitively 
on the small parameters uj p /ujb and 9. As a guide we can use a much simpler 
condition for anomalous cyclotron resonance on X mode, keeping in mind that X 
mode has phase velocities larger than L-0 mode. In the limit w<Cwg resonance 
on X mode occurs at frequencies 

u~J# (8) 

Note, that for 7 ^> (ujb/^ p ) 2 both the resonant frequency and frequency dif- 
ferences are much smaller than cyclotron frequency. In case of relativistically 
streaming plasma on open field lines, the rhs of Eq. (|8j) is multiplied by 7 3 
(where 7 P is Lorentz factor of plasma bulk motion), as plasma density, wave 
frequency and Lorentz factor are all smaller by ~ 7 P . As a result, for the fast, 
high field pulsar like Crab the resonance condition for radio waves cannot be 
satisfied inside the magnetosphere for reasonable plasma parameters. From (JsJ) 
it follows that in order to have resonant frequency below ~ 10 GHz, it is required 
that 7A > 10 13 . This provides an order-of-magnitude estimate for the required 
plasma parameters. In making this estimate we neglected the Lorentz factor 
arising due to rotational velocities of emitting plasma, which can be important 
near the light cylinder. 

Using the value of the magnetic field at the light cylinder and parametriza- 
tion of plasma density ^ we search through three parameters: plasma density 
normalization A, Lorentz factor of the fast particles 7 and viewing angle 9. To 



fit the results, we use Fig. 4 of Eilek & Hankins (20071, trying to reproduce 
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the frequencies and separations of the emission bands. We find a satisfying fit 
for the following parameters, see Fig. [4] Plasma over-density A = 3.3 x 10 5 
(note, that this value is close to what was inferred for the plasma density on 
closed field lines of pulsar B in the double pulsar PSR J0737— 3039 (Lyutikov 
& Thompson 2005 1). Lorentz factor of the fast beam is 7 = 7.4 x 10 Y ; this 



turns out to be close to the estimate of radiation- limited acceleration (10 1. The 
viewing angle with respect to magnetic field is 8 = 0.0022. A reasonable fit is 
achieved in a fairly limited volume of parameter space, 3.2 x 10 5 < A < 3.8 x 10 5 , 
6.9 x 10 7 < 7 < 7.4 x 10 7 , 0.00251 < < 0.00225. Most importantly, for the 
first four bands the value of bands' separation 2(v s+1 — v s )/{ v s+i + v s ) are 
.036, 0.053, 0.077, which is very close to the mean quoted value of 0.06. 



3.2 Growth rates 

Since for our best fit parameters 9 ^> (u) p /ujb) 2 , Eq. ^ implies the L-0 mode 
is nearly transverse at resonance. We can then use the previously calculated 
growth rates for excitation of transverse electromagnetic modes at anomalous 



cyclotron resonance (Kazbegi et al. 1991 Lyutikov et al. 19991 



p,res 



(9) 



where co Pires is the plasma frequency of the fast beam, generated during recon- 
nection. Note, that since the emission bands' spacing is much smaller than their 
carrier frequencies, the growth rates are nearly the same for different bands (ob- 
served intensity also depends on saturation mechanism). Similar growth rate 
and, presumably, similar saturation levels explain generation of multiple emis- 
sion bands. We expect that during a reconnection event, all the particles in the 



reconnection region may be accelerated, 
rate can be very high r ~ lj p 



so that the Growth 



3.3 Lorentz factor of fast particles on closed lines 

We envision that reconnection events close to the Y-point (near the last closed 
field lines at magnetic equator) accelerate particles along the closed field lines. 
The maximum available potential is related to the resistive process occurring 
during field reconnection, like the speed of reconnection and a number of flux 
tubes being reconnected. In addition, for nearly orthogonal rotators the "null 



charge surfaces" Cheng & Ruderman (19771, where Goldreich-Julian density 



vanishes, lies inside the polar cap, so that after reconnection a large potential 
drop would develop along the field lines Istomin ( 2004 ) . The upper bound, which 



is not reached by far, can be estimated as the total potential across the open field 
lines, $ ~ B N sRns(QR/ c ) 2 ', this corresponds to a Lorentz factor j max ~ 10 11 
( Istomin| ( |2004[ ) esimated Lorentz factor is about an order of magnitude smaller). 
The maximum electron energy will, in fact, be limited by radiative damping, like 
curvature and IC radiation. Assuming that curvature radiation is the dominant 
loss process, we can estimate the maximum energy than an electron can reach 
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in accelerating electric field E ~ B (assuming that reconnection inflow velocity 
is close to the speed of light). We find 

/B\ 1/4 /— / r \~ 1/4 

7~( 7 J ^-6.7X10^— ) (10) 

where R c ~ c/f2 is the curvature radius. Surprisingly, this matches nearly ex- 
actly the value that we inferred from fitting the observed narrow-band structure. 



3.4 High plasma density on closed field lines 



According to the model, a high density plasma should be present on closed field 
lines. The overdensity we find here, modeling Crab GPs, ~ 10 5 , turns out to 
be close to the one found by eclipse modeling of the binary pulsar system PSR 
J0737-3039A/B Lyutikov & Thompson (20051. How related are these peices 
of evidence and how dense plasma is created and supported on closed field 
lines? These two cases seem to very different, in fact. In slow pulsars, like 
PSR J0737-3039B, plasma can be efficiently trapped on closed field lines by 
magnetic bottling (for ~ 10 6 periods in case of PSR J0737-3039B, Lyutikov & 



Thompson (2005J), but not in case of Crab, in w hich case radiative decay times 



are too short, ~ 10 sec at the light cylinder. Lyutikov & Thompson (20051 



proposed that high densities on closed field lines of PSR J0737-3039B may be 
explained by interaction with the wind of the companion, but similar overdensity 
in an isolated pulsar questions that possibility. Though there are several ways 
particles can populate closed field lines {e.g. kinetic drift from open field lines, 
pair production by the 'backward" beam from outer gaps), the demands in case 
of short pariod pulsars like Crab are pretty high: the mechanism should create 
an over-density of the order 10 5 with no efficient bottling. 

In case of Crab this over-density is required in a fairly limited region near the 
light cyllinder, where IGPs are presumably produced. It is somewhat natural 
to associated this density enhancement with the Y-point, where the last closed 
field line approaches the magnetic equator. In fact, we indeed may expect high 
density around that region. First, even in rigidly rotating dipolar magnetosphere 
the charge density diverges at that point [Goldreich fc Julian (19691. Second, 
in case of a 



more realistic force-free aligned rotator, poloidal magnetic field 
also diverges at the Y-point Gruzinov (2005); Spitkovsky (20061. To resolve 



this divergency, non-ideal effects such as resistivity Spitkovsky (20061 and/or 
particle inertia Komissarov ( 2006 1 should be taken into account. At the moment, 



we leave this possibility to further studies. 



3.5 Predictions: simultaneous GLAST signal 

The model has interesting prediction for the forthcoming GLAST mission. The 
high energy beam is expected to produce curvature radiation at ~ hj 3 £l ~ 
30 GeV. This energy corresponds to the maximum effective area of the LAT 
instrument. Thus, one expects to see GeV photons contemporaneous with GPs. 
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Estimating the beam density ~ uqj and emitting volume as ~ 0.1R^ C , the 
total power in curvature photons is ~ 2 x 10 3 6 erg/sec, or 4 x 10 3 7 photons/sec. 
Assuming isotropic emission (this gives a lower limit on the observed flux), the 
expected flux at Earth is ~ 10~ 7 photons cm~ 2 s _1 during a GP. This is nearly 
three orders of magnitude higher than the integral sensitivity of LAT at 30 GeV 
(for high latitude sources; Galactic background will also contribute in case of 
Crab), implying that a duty cycle of > 10~ 3 for GP emission is needed to detect 
GP in GeV range. We conclude that we may reasonably expect GeV photons 
contemporaneous with GPs. 



4 Implications of the model 

Recent advances in observational technique uncovered very detailed properties 



of pulsar radio emission, and of Crab giant pulses in particular (Eilek & Hankins 



20071. Inspired by this work, in the present paper we describe a quantitative 
model of Crab giant pulses. To the best of our knowledge this is the first work 
on pulsar emission mechanism that actually does a fitting of pulsar spectra 
to a particular model. The model can reproduce fairly well the narrowband, 
unequally spaced emission bands seen in GHz frequency range in Crab giant 
pulses associated with the interpulse. It also explains why banded structure in 
IGPs in Crab is not seen below ~ 4 GHz: anomalous cyclotron resonance on 
L-0 mode occurs above plasma frequency, which for our typical parameters is 
in the GHz range. Pulse to pulse variation in location of emission bands are 
due to fluctuations of the plasma density. 

Note, that though the emission involves cyclotron transitions, both the emis- 
sion frequencies and inter-band spacing are much smaller than the local cy- 
clotron frequency for the best fit parameters. This is a curious property of 
anomalous cyclotron resonance. 

Perhaps the most striking features of the model is that GPs are generated on 
closed field lines, in stark contrast to all other models of pulsar radio emission. 
We stress that this only applies to GPs, regular pulses are generated on open field 
lines, as is well established by a multitude of observational facts. We argue that 
giant pulses are different. The fact that we see GPs in some particular pulsars 
and that only the Crab IGP shows a banded frequency structure, is, in some 
sense, just a chance occurrence. The resonance condition is a sensitive function 
of the angle between the line of sight, local magnetic field and plasma density, 
so that a mismatch by ~ 10~ 3 radians in observer angle cither destroys the 
resonance completely or places the resonant frequency outside of the observed 
frequencies. This explains why GPs are typically seen in narrow phase windows. 
This may also be related to the fact that only few pulsars show GP phenomenon: 
in those cases we are lucky in terms of orientation, so that the line of sight 
passes close to the direction of local magnetic field at the reconnection region, 
and in terms of local plasma parameters, so that the emission frequency, which 
scales nearly cx ujg falls in the observed frequency range. Exact location of the 
reconnection region is not specified in the model, but it is naturally to assume 
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that it occurs close to the Y point, where we expect high current concentrations 
( |Gruzinov[[2005| . 

One of the main implications of the model is that (at least near the Y-point) 
closed field lines are populated with dense plasma, exceeding by a factor 10 5 — 10 6 
the minimum charge density. On the one hand, this is not expected from the 
basic pulsar model (Goldreich & Julian 19691. On the other hand, presence 



of such dense plasma on closed field lines is solidly established in the case of 
the double pulsar PSR J0737-3039 ( [Lyutikov fc Thompson} |2005| ). Somewhat 
surprisingly, the over-density ~ 10 5 (with respect to Goldreich- Julian density) 
we find here, modeling the generation of Crab giant pulsars, is similar to the 
over-density required to explain eclipses in the double pulsar. 

To fit our findings into a global model of pulsar radio emission, we note 
that at frequencies above 2 GHz (but below NIR), Crab emission around the 
interpulse undergoes drastic changes with respect to both lower radio frequen- 
cies and higher optical through X-ray emission Moffett & Hankins (19961: the 
interpulse is shifted by ~ 10 degrees and new High Frequency Components 
(HFC1 and HFC2) trailing the interpulse appear. This is the range where the 
banded structure appears in IGPs. Our findings are consistent with the pos- 
sibility that interpulse emission at these frequencies is a qualitatively different 
component from the main pulse and from the low frequency interpulse. Still a 
clear (geometrical?) picture of emission components is still missing. 

For small angles of propagation with respect to magnetic fields , the disper- 
sions of the X and the L-0 modes nearly coincide, while polarization of both 
are nearly orthogonal for the parameters of interest. We have shown that the 
pattern of emission bands seen in IGPs matches well the anomalous cyclotron 
resonance on the L-0 mode. On the other hand, resonance bands of the X 
mode may lie fairly close to the L-0 modes, so it is natural to expect that the 
X mode should also be emitted. Emission bands associated with the X modes 
have distinctly different properties from L-0 mode, being spaced in nearly equal 
intervals, equal to the frequency of the first harmonic. Large spacing, nearly 
the same as carrier frequency, makes it nearly impossible to observe X-mode 
emission bands with a finite receiver bandwidth (additional complications at 
low frequncies come from interstellar scattering) . Still, the fact that a spectrum 
is not continuous, but consisting of emission bands, may be verified statistically 
during simultaneous observations of GPs in different frequency bands. Such ob- 
servations of the millisecond pulsar B1937+21 (Popov & Stappers 2003 1 showed 
that out of 10-15 GPs observed at two frequency windows 14141446 MHz and 
22102250 MHz, no events were found to occur simultaneously at both frequen- 
cies. This requires that spectra of individual GPs are narrow band, /S.v/v < 0.5, 
consistent with our picture. 

Additional complications may come from absorption of the X-mode at cy- 
clotron resonance. Since at the light cylinder the radiation frequency is much 
smaller than the local cyclotron frequency, resonant absorption may occur only 
in the wind (Sincell & Krolik 19921. At distances much larger than the light 



cylinder, the radially propagating L-0 mode is polarized nearly along magnetic 
field and thus is not absorbed at the cyclotron resonance ( Petrova 2006 1 , while 
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the X mode may experience absorption. 

The present model can also accommodate qualitatively the fact that some- 
times emission bands seem to be drifting slightly up in frequency. If the fast 
beam is generated close to the magnetic equator, where magnetic field is the 
weakest, then as it propagates along the closed field lines, the local cyclotron 
frequency increases, leading to an increase of the observed frequency. Thus, plac- 
ing beam generation region at magnetic equator explains why emission bands 
always drift up in frequency (Eilek & Hankins, priv. comm.). 

At the moment it is not clear how rare or ubiquitous is the phenomenon 
of GPs in pulsars. As many as tens percent may show GPs (Ransom, priv. 
comm.). Of those pulsars that do show GPs, detection of banded structure in 
GPs depends sensitively on the plasma parameters close to the light cylinder 
and on the angle between line of sight and local magnetic field . Taking our best 
fit value of 6 ~ 2 x 10 -3 we can estimate a relative number of pulsars in which 
we expect to see bands in GPs as ~ 9. Such an estimate is very rough and is 
naturally subject to a number of uncertainties, like beaming fraction of normal 
and giant pulses. 

Finally, we note that closed field lines of the Earth magnetosphere are very 
active in producing high brightness radio emission like auroral hiss, roars and 



burst (LaBelle & Treumann 2002), though they all occur at normal cyclotron 
resonance (and thus require a loss cone-type distribution) at frequencies of the 
order of the local cyclotron frequency (and thus are too high to be directly 
applicable to observations of Crab pulsar at radio waves). It is intriguing to 
speculate that similar type emission might be observed from Crab at ~ 10 12 — 
10 13 Hz, but the required observations at such frequencies are challenging. 

I would like to thank Timothy Hankins, Jean Eilek, Qinghuan Luo, George 
Machabeli, Alison Mansheim, Mikhail Popov and Scott Ransom for discussions 
that were vital for this work. 
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A Other excitation mechanisms 



In this appendix we discuss other excitation mechanisms and point out why 
they cannot be responsible for Crab GPs, at least for those GPs that show a 
narrowband structure. In doing so, we assume that the observed band structure 
is due to emission processes, not propagation. 

At uj <C ujb there are three modes: Alfven , X and L-O, which can be excited 
at Cherenkov (except X mode), cyclotron and anomalous cyclotron resonances. 
We have discussed excitation of the X and the L-0 modes at anomalous cy- 
clotron resonances; let us now discuss other possibilities. 



A.l Normal cyclotron resonance 

For instability to operate on normal cyclotron resonance particle distribution 
should be of the loss-cone type. Such distributions may be expected in the 
outer gap models, which produce secondary particles with non-zero transverse 
momenta. Effects of magnetic bottling on a down-ward propagating particles 
will create a loss-cone. In addition, since an outer gap operates near the last 
closed field lines, it is expected that some particles will drift from open to 
closed lines. This possibility is attractive since the closed field lines may store 
particles bouncing between the poles for a long time. Development of the loss- 
cone instability would then convert this stored energy into radiation. This would 
create radio emission from closed field lines. Both the X and the L-0 modes 
can be emitted at cyclotron resonance, but the main problem is that emitted 
frequencies are typically close to cyclotron frequency ujb, which is some three 
orders of magnitude higher than observed. In principle, the resonant frequency 
can be smaller than lob as, e.g. , in case of backward wave oscillators, but this 
would require fine tuning the frequencies in one part in a thousand. In addition, 
inter-band frequency spacing is large, ~ub. 

A. 2 AlfVen mode 



Excitation of Alfven waves in pulsars has been considered by Lyutikov ( 2000 1 



Alfven waves exist for < ui < ^| ( ui B + 2uj p — y (wfj + 2uj p ) 2 — 8u B ujp cos 2 9 

^/2uj p cos9. At low frequencies, u> <C uj p cos9, uj p <C ujb, dispersion relation be- 
comes 

/ uj 2 (kc) 2 sin 2 o\ 
uj = kccos9 1 1- ^ , for uj <C uj p cos9, uj p <C lob (11) 

Alfven waves can be excited at Cherenkov resonance at 



] 



sm9 y % 0J% 
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which requires ojb / '(V^^pjb) > !• The resonant frequency is 



UJ 



/2uj p 

76 



cot / 



(13) 



so that at a given angle emission is generated with given uj. Though the reso- 
nance condition can, in general, be satisfied, it not clear how to produce different 
emission bands. 

Alfven waves can also be emitted at anomalous Doppler resonance. For small 
angles, 9 <C ^ p kc , the resonant frequencies are similar to those of X mode, while 



for 9 > 



we find 



SUJbUJ 



V 



7b cos 9 sin 9 



1/3 



(14) 



In addition, condition uj < u p requires 



> 



(15) 



An appealing property of the resonant condition) 14 1 is that it produces a set 



of frequency bands, separated by frequencies much smaller than cyclotron fre- 
quency lob- On the other hand, inter-band separation decreases with fre- 
quency, see Fig [5] contrary to observations. This invalidates the model of |Is-| 
(20041, at least in application to IGPs. In addition, Alfven waves cannot 



tomin 



escape from plasma and need to be converted into escaping modes before they 
are damped on plasma particles. 
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Figure 2: Wave dispersions u(k) in pair plasma in strong magnetic field, 
u>b 3> w p , for oblique propagation. At low frequencies u> <C ujb there are 
three modes labeled X (polarized orthogonally to k -B plane), Alfven and 
L-0 (both polarized in the k -B plane). The L-0 mode has a resonance at 
~ u>b + sin 2 8/ujb and cut-off at \/2uj p . The Alfven mode has a resonance 
at ~ v^^pCos^. The sign || indicates locations where corresponding waves 
are nearly longitudinally polarized. The two high frequency, uj > ujb, waves 
with nearly identical dispersion have a cut-off at the upper hybrid frequency 

ujjh = \/ lj2 b + 2w 2 . 
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Figure 3: Resonant frequencies for the L-0 (solid line) and the X modes (doted 
line) for anomalous cyclotron resonance as function of angle between the line of 
sight and magnetic field for first four resonances s = —1, ... — 4. The cut-off of 
the L-0 mode at 8 = is at total plasma frequency \f2io v j (2-n) = 7.76 GHz. 

The X modes bands are equally spaced at 9 — with 5v — (corresponding 

to 1.52 GHz); higher order resonances are not shown. The vertical dashed line 
gives the best fit values of 8 — 0.0022 for resonances on L-0 mode. 
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Figure 4: Location of emission bands (white stripes) for the fiducial model. The 
fitted observations correspond to Fig. 4 of Eilek & Hankins (20071. 
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Figure 5: Resonant frequencies for anomalous excitation of Alfven waves. Pa- 
rameters are the same as in Fig. [4j Alfven waves frequencies are limited to 
oj < y/2ui p cos 9. Contrary to observations, the inter-band spacing decreases 
with frequency. 
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